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Abstract The present work deals with determining the 
necessary parameters considering a three dimensional 
model to simulate in a realistic way the turning process 
on machine tool. 

This paper is dedicated to the study of the self-excited 
vibrations incidence on various major mechanics charac- 
teristics of the system workpiece / tool / material. The 
efforts (forces and moments) measurement using a six 
components dynamometer confirms the tool tip moments 
existence. 

The fundamental frequency of 190 Hz proves to be 
common to the tool tip point displacements, the action 
application point or at the torque exerted to the tool tip 
point. The confrontation of the results concerning dis- 
placements and efforts shows that the applications points 
of these elements evolve according to similar ellipses lo- 
cated in quasi identical planes. The large and the small 
axes of these ellipses are increasing with the feed rate 
motion values accordingly to the mechanical power in- 
jected into the system. Conversely, the respective axes 
ratios of these ellipses are decreasing functions of the 
feed rate while the ratio of these ratios remains constant 
when the feed rate value is increasing. 

In addition, some chip characteristics are given, like 
the thickness variations, the width or the hardening phe- 
nomenon. 
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Nomenclature 



A 

[A]o 
ap 

a f ( b f) 



BT 
BW 



f 

fcop 

F v (F n ) 
F x 

Fy 

F z 

hmax 

h 

M A 
M Q 
N 

71/ 

n fa (n fb ) 

o 

O' 
Pf 



Central axis point; 

Actions torsor exerted on the tool tip in 
O point; 

Cutting depth (mm); 
Large (small) ellipse axis attached to 
the geometrical place of forces applica- 
tion points; 
Block Tool; 
Block Workpiece; 

Directory line of moments projection to 

the central axis (i = 1, 2); 

Ellipse point belong attached of forces 

application points (i=l-5); 

Feed rate (mm/rev); 

Chip segmentation frequency (Hz); 

Variable (nominal) cutting force (N); 

Effort on cross direction (N); 

Effort on cutting axis (N); 

Effort on feed rate axis (N); 

Maximum chip thickness (mm); 

Minimum chip thickness (mm); 

Chip undulation length (mm); 

Cutting forces minimum moment in A 

who actuate on the tool (N.m); 

Cutting forces moment in the O point 

who actuate on the tool (N.m); 

Spindle speed (rpm); 

Pf normal direction; 

Normal direction projection to the Pf 

on the af (bf); 

Tool point reference; 

Dynamometer center transducer; 

Plane attached to the forces application 

points; 
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Nomenclature continuation 



P u Plane attached to the tool point displace- 

ments; 

R Cutting forces vector sum who actuate on the 

tool (N); 

r e Cutting edge radius (mm); 

R Sharpness radius (mm); 

R t Roughness (/xra); 

T Time related to one revolution of workpiece 

(s); 

u Tool tip point displacement (m); 

V Cutting speed (m/min); 

w max Maximum chip width (mm); 

w m i n Minimum chip width (mm); 

WTM Workpiece-Tool-Machine tool; 

x Radial direction; 

y Cutting axis; 

z Feed rate direction; 

a Clearance angle (degree); 

a K(xy) Angle of main tool displacements direction 
included in the plane (x,y) (degree); 

C£ K ( yz ) Angle of main tool displacements direction 
included in the plane (y,z) (degree); 

At Time corresponding of phase difference be- 

tween two signals (s); 

<P Primary share angle (degree); 

(f c Chip width slopes angle between each undu- 

lation (degree); 

iff Ui Phase difference between the tool tip dis- 
placements components i=x, y, z (degree); 

7 Cutting angle (degree); 

A s Inclination angle of tool edge (degree); 

K r Direct angle (degree); 

@e(xy) Stiffness principal direction angle related to 

the plane (x,y) (degree); 
@e(yz) Stiffness principal direction angle related to 

the plane (y,z) (degree); 
£ c Chip hardening coefficient; 



1 Introduction 

In the three-dimensional cutting case, the mechanical ac- 
tions torsor (forces and moments), is often truncated: the 
moments part of this torsor is neglected fault of adapted 
metrology [27]. However, efforts and pure moments (or 
torque) can be measured [12J. Recently, an application 
consisting in six components measurements of the ac- 
tions torsor in cutting process was carried out for the case 
of high speed milling [14], drilling [23], [35], etc. Cahuc 
et al., in [9] presents another use of this six components 
dynamometer in an experimental study: the taking into 
account of the cut moments allows a better machine tool 
power consumption evaluation. The present paper is ded- 
icated to the six components dynamometer use to reach 
with fine accuracy a vibrations dynamic influence evalu- 



ation in turning on the system Workpiece/ Tool/ Machine 
tool (WTM). 

The concepts of slide block system, moments and 
torsor (forces and moments) are directly related to the 
work, undertaken to the previous century beginning, on 
the mathematical tool "torsor" [29]. Unfortunately, until 
now the results on the cutting forces are almost still vali- 
dated using platforms of forces (dynamometers) measur- 
ing the three components of those [24J, [26J. The actions 
torsor is thus often truncated because the torsor moment 
part is probably neglected fault of access to an adapted 
metrology [33]. [34]. 

However, recent experimental studies showed the ex- 
istence of cutting moments to the tool tip [10], [21J, [22J, 
[32] . Their taking into account allows thus the best ma- 
chine tools output approaches [9], [15]. Nowadays, the 
use of a dynamometer measuring the mechanical actions 
torsor six components [4], [13], [23] allows a better cut 
approach and should enable to reach new system WTM 
vibrations properties in the dynamic case. 

Moreover, the tool torsor has the advantage of being 
transportable in any space point and in particular at 
the tool tip in O point. The study which follows about 
the cut torsor carries out in several stages including two 
major; the first relates to the forces analysis and second 
is dedicated to a first moments analysis to the tool tip 
during the cut. 

After this general information on the torsor concept, 
an analysis of the efforts (section [2]) exerted in the cut- 
ting action is carried out. Thanks to this analysis and 
experimental tests, we establish that the points of vari- 
able application exerted on the tool type are included in a 
plane, more precisely following an ellipse. Then, accord- 
ingly in section [3j the first moments analysis obtained at 
the tool tip in O point is performed. The torsor central 
axis is required (section [4)) and the central axes beams 
deduced from the multiple tests strongly confirm the mo- 
ments presence to the tool tip. 

Some chip characteristics are presented in the sec- 
tional and compaired to experimental data base. Before 
concluding, a correlation (section EJ) displacements / cut- 
ting forces shows that these two dual elements evolve in 
a similar way following ellipses having the similar prop- 
erties. 



2 Forces analysis 

2.1 Tests results 

The experiments are performed within a framework sim- 
ilar to that of Cahuc et al [10J. For each test, except the 
feed rate values, all the turning parameters are constant. 
The mechanical actions are measured according to the 
feed rate (f) using the six components dynamometer [12] 
following the method initiated in Toulouse [31], devel- 
oped and finalized by Couetard [13] and used in several 
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occasions [9], [14], [15], [21], [22]. On the experimental 
device (Fig. [I]) the instantaneous spindle speed is perma- 
nently controlled (with an accuracy of 1%) by a rotary 
encoder directly coupled with the workpiece. During the 
tests the insert tool used is type TNMA 16 04 12 carbide 
not covered, without chip breeze. The machined mate- 
rial is an alloy of chrome molybdenum type 42CrMo24. 
The test-workpieces are cylindrical with a diameter of 
120 mm and a length of 30 mm. They were designed 
starting from the Finite Elements Method being coupled 
to a procedure of optimization described in [6]. 




Fig. 1 Experimental device and associated measurement el- 
ements. 



Moreover, the tool geometry is characterized by the 
cutting angle 7, the clearance angle a, the inclination 
angle of edge j s , the direct angle n r , the cutting edge 
radius r e , and the sharpness radius R [20]. In order to 
limit to the wear appearance maximum along the cutting 
face, the tool insert is examined after each test and is 
changed if necessary (Vb < 0.2 mm ISO 3685). The tool 
parameters are detailed in the Table [H 



the lowest value is the (F x ) one, while the highest in the 
absolute value is (the F z ) one. While taking as reference 
the absolute value of (F x ) the following relation between 
these three components comes: \F X \ = = 

In the unstable case we observe that the efforts com- 
ponents on the cutting axis (F y ) has the most important 
average amplitude (1 500 N). It is also the most disturbed 
(± 700 N) with oscillations between -2 200 N and -800 N. 
In the same way the effort according to the feed rate axis 
(F z ) has important average amplitude (1 000 N) and the 
oscillations have less width in absolute value (± 200 N) 
but also in relative value (± 20%). As for the effort on the 
radial direction (F x ) it is weakest on average (200 N) but 
also most disturbed in relative value (± 200 N). These 
important oscillations are the tangible consequence of 
the contact tool/ workpiece frequent ruptures and thus 
demonstrate the vibration and dynamical behaviour of 
the system WTM. 

Finally, we note that the amplitudes of all these ef- 
forts components applied to the tool tip are slightly de- 
creasing functions of time in particular for the compo- 
nent according to the cutting axis. 
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Fig. 2 Signals related to the resultant components of cutting 
forces following the three cutting directions; test case using 
parameters ap = 2 mm, f = 0,1 mm/rev and N = 690 rpm. 
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a 






r e 


R 


-6" 


6" 


-6" 


91" 


1,2 mm 


0,02 mm 



Table 1 Tool geometrical characteristics. 



Two examples of resultant efforts measurements ap- 
plied to the tool tip are presented: one of these for the 
stable case, ap = 2 mm (Fig. [2]), and other for the case 
with instability, ap = 5 mm (Fig. [3]). In the stable case 
it appears that the force components amplitudes remain 
almost independant from time parameter. Thus, the am- 
plitude variation is limited to 1 or 2 N around their nom- 
inal values, starting with 200 N for (F x ) and until 600 N 
for (F y ). These variations are quite negligible. Indeed 
the nominal stress reached, the component noticed as 




Fig. 3 Signals related to the resultant components of cutting 
forces following the three (x, y, z ) cutting directions; test 
case using parameters ap = 5 mm, f = 0,1 mm/rev and N = 
690 rpm. 
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2.2 Frequency analysis 

The signals frequency analysis performed by using FFT 
function enables to note in Fig.[4]the presence of frequen- 
cies peaks around 190 Hz. Around this frequency peak, 
we note for the three forces components, a quite high 
concentration of energy in a wide bandwidth around 70 
Hz (36% of the fundamental frequency). All things con- 
sidered, this width of frequency is of the same order of 
magnitude as observed (13% of fundamental) by Dimla 
[T6] for a depth of cut three times lower (ap = 1.5 mm) 
but for an identical feed rate (f=0.1 mm/rev) and a cut- 
ting speed similar. This remark confirms that the efforts 
components is proportionnal to the depth of cut ap as 
indicated by Benardos et al., [I]. 
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Fig. 5 Cutting force F v evolution around the nominal value 
F n ; test case using cutting parameters ap = 5 mm, f = 
0,1 mm/rev and N = 690 rpm. 
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Fig. 6 The resultant cutting force components variable zoom 
on the three cutting directions; test case using parameters ap 
= 5 mm, f = 0,1 mm/rev and N = 690 rpm. 



Fig. 4 Cutting forces magnitude FFT signal on the three 
(x, y, z) directions; test case using cutting parameters ap = 
5 mm, f = 0,1 mm/rev and N = 690 rpm. 



Let us recall that the same frequency of 190 Hz was 
observed in the tool tip displacements case (in confor- 
mity with [4]). Consequently, the cutting forces compo- 
nents variations and the self-excited vibrations are influ- 
enced mutually, in agreement with [17], [IS], [25J, |30j. 
Also, in agreement with research on the dynamic cut- 
ting process [28], we note that the self-excited vibrations 
frequency is different from the workpiece rotational fre- 
quency which is located around 220 Hz. 



2.3 Forces decomposition 

The forces resultant components detailed analysis high- 
lights a plane in which evolves a variable cutting force F v 
around a nominal value F n (see further). This variable 
force is an oscillating action (Fig. [5j) which generates u 
tool tip displacements and maintains the vibrations of 
elastic system block-tool BT [7]. 

Thus, the cutting force variable (Fig. [5]) and the self- 
excited vibrations of elastic system WTM are interac- 
tive, in agreement with research work [18J, |28j. 



The cutting forces variable part can be observed and 
compared. Not to weigh down this part, the cutting forces 
analysis is voluntarily below restricted at only two dif- 
ferent situations: 

- stable process using cutting depth ap = 2 mm (Fig.[7fei), 

- unstable process (with vibrations) using cutting depth 
ap = 5 mm (Fig. 0)) 



i\p = 2 mm 
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Fig. 7 Stable process (a) and unstable process (b). 



Self-excited vibrations in turning : Forces torsor analysis 



5 



The vibrations effects on the variable forces evolution 
are detailed on the Fig. [3d. Moreover, the variable forces 
and displacements analysis associated with the tool tip 
at the time of the unstable process shows that the forces 
variation ratio is equivalent to the tool tip displacements 
variation ratio (in conformity with [2j). This aspect will 
be quantified further. Also we concentrate our attention 
mainly on the unstable case (ap = 5 mm). 



2.4 Plane determination attached of the forces 
application points 

The tests analysis shows that, in the vibratory mode, the 
load application points describe an ellipse (Fig. [8]), that 
is not the case in the stable mode (without vibrations, 
Fig. Ei)). 

The method used in [3] to determine the tool tip dis- 
placements plane is taken again here (in conformity with 
Appendix, section l8T|) to establish the plane Pf, place of 
the load application points, characterized by its normal 
n f (Table El). 



, 1000 








•HI 



Fig. 8 The ellipse plane Pf attached of range of the forces 
application points considering ap = 5 mm, f = 0,1 mm/rev 
and N = 690 rpm. 



ap = 5 mm 


X 


y 


z 


Normal 


n f 




0.1 


0.46 


- 0.1 


- 0.882 


f (mm/rev) 


0.075 


0.419 


- 0.097 


- 0.903 


0.0625 


0.292 


- 0.113 


- 0.95 




0.05 


0.245 


- 0.107 


- 0.964 



Table 2 Directory normal of the plane Pf attached of forces 
application points on the tool; case study using ap = 5 mm, 
f = 0,1 mm/rev and N = 690 rpm. 



As for the tool tip displacements study related, a new 
reference system (nf a ,nfb) is associated at the load ap- 
plication points ellipse place. In this new reference sys- 
tem (in conformity with section [8^2]) the ellipse large (re- 
spectively small) axis dimensions a/ (respectively bf) are 
obtained (Fig. [9]). The values of a/, bf as that of their 
ratio df I bf are consigned in the Table 03 It should be 
noted that cif and bf are increasing feed rate (/) func- 
tions. It is thus the same for the ellipse surface which 
grows with (/), in perfect agreement with the mechanical 
power injected into the system which is also an increas- 
ing feed rate function. On the other hand, the ratio a/ 
/ bf is a decreasing feed rate function. Thus the ellipse 
elongation evolves in a coherent way in dependence with 
the feed rate. 





— 1500 




— 900 




-300 £ 
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— -30Q 





Fig. 9 Ellipse approximation attached of forces applica- 
tion points; case study using parameters ap = 5 mm, f = 
0,1 mm/rev and N = 690 rpm. 



/ (mm/rev) 


a f (N) 


b f (N) 


af/bf 


0.1 


690 


280 


2.46 


0.075 


545 


220 


2.47 


0.0625 


365 


147 


2.48 


0.05 


220 


88 


2.5 



Table 3 Large and small ellipse axes attached of forces ap- 
plication points depending on feed rate parameter; case study 
using ap = 5 mm, f = 0.1 mm/rev and N = 690 rpm. 



Let us look at now the moments evolution to the tool 

tip. 
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3 First moment analysis 

3.1 Experimental results 

For each test, the mechanical actions complete torsor 
is measured according to the method already detailed 
in the section I27T1 [T3] . Measurements are taken in the 
six components dynamometer transducer O' center and 
then transported to the tool point O via the moment 
transport traditional relations |8j. As for the forces, the 
moments variable part is extracted from measurements. 
An example of the results measurement is given on the 
Fig. [10] which zooms on the moments variable part is 
presented in the Fig. [Til 



Znom (Figure 1 1 ) 




Fig. 10 The moments components time signals that action 
on the tool tip following the three directions; case study using 
ap = 5 mm, f = 0.1 mm/rev and N = 690 rpm. 



Taking into account the recordings chaotic aspect 
obtained, an accurate moments components frequency 
analysis is necessary. It is the object of the following sec- 
tion. 



3.2 Moments frequency analysis 

An example of moments signals frequency analysis dur- 
ing the vibratory cutting is presented in the Fig. [T2j As 
for the forces analysis, the moments components FFT 
shows that the most important frequency peak is local- 
ized around 190 Hz. 

Moreover of all the components, the most important 
fundamental amplitude is that corresponding to the mo- 
ments components following to the cutting axis (y) as its 




Fig. 11 The moments components variable part zoom that 
action on the tool tip following the three cutting directions; 
case study ap = 5 mm, f = 0,1 mm/rev and N = 690 rpm. 



transport at the tool tip confirms it. It should be noted 
that the force component following to this same axis is 
also most important but has obviously no influence on 
the moments related to this axis due to co-linearity of 
these two elements. Conversely, the least important vi- 
bration amplitude is that of the moment component lo- 
cated on z axis. However, the number of revolutions be- 
ing important according to this axis the M z component 
contribution to the torques power consumption remains 
significant. 




Fig. 12 Signal FFT related to the moment on the three (x, 
y, z) directions; case study using ap = 5 mm, f = 0.1 mm/rev 
and N = 690 rpm. 



The appearance of other peaks, which are harmonics, 
slightly modifies the three-dimensional moments repre- 
sentation (Fig. [13]) which is not exclusively any more in 
a plane although the ellipse essence is in a plane. This 
representation approaches a light form of letter eight con- 
trary to the elliptic planar form efforts characteristic. 
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,-200 




Fig. 13 Moments place space representation. 

4 Central axis 

4.1 Central axis determination 

It is well-known that with any torsor, it is possible to 
associate a central axis (except the torsor of pure mo- 
ment), which is the single object calculated starting from 
the torsor six components [8]. 

A torsor [A] Q in a point O is composed of a forces 
resultant R and the resulting moment Mq : 



[A] 



f R, 



0_ \Mo. 

The central axis is the line defined classically by: 



(i) 



. R A M , ,, 



(2) 



where O is the point where the mechanical actions 
torsor was moved (here the tool tip) and A the current 
point describing the central axis. OA is thus the vector 
associated with the bi-point [O, A] (Fig. [T4|) . 

This line (Fig. [T4l-(a)) corresponds to points geomet- 
ric place where the mechanical actions moment torsor 
is minimal. The central axis calculation consists in de- 
termining the points assembly (a line) where the torsor 
can be expressed according to a slide block (straight line 
direction) and the pure moment (or torque) [8]. 

The central axis is also the points place where the re- 
sultant cutting force is co-linear with the minimum me- 
chanical moment (pure torque). The test results enable 




Fig. 14 Central axis representation (a) and of the co- 
linearity between vector sum R and minimum moment Ma 
on central axis (b). 



to check for each point of measurement the co-linearity 
between the resultant cutting force R and moment Ma 
calculated related to the central axis (Fig. [14]- (b)). 

The meticulous examination of the mechanical ac- 
tions torsor six components shows that the forces and 
the moment average values are not null. For each mea- 
sure point, the central axis is calculated, in the stable 
(Fig. [T5l-(a)) and unstable mode (Fig. H5]-(b)). In any 
rigour the case ap = 2 mm should be described as quasi- 
stable movement, because the vibrations exist but their 
amplitudes are very low - of the order of the fim -, thus 
quasi null compared to the other studied cases. Consid- 
ering the cutting depth value ap = 3 mm, the recorded 
amplitude was 10 times more important. 



up = 2 

f= 0.0625 inin rev 



sip = 5 inin 
f = 0.0625 mm/rev 




Fig. 15 Central axes representation obtained for 68 rpm the 
workpiece speed and feed rate f = 0,0625 mm/rev; a) - stable 
process ap = 2 mm; b) - unstable process ap = 5 mm. 



In the presence of vibrations (ap = 5 mm) for a 
68 rpm the workpiece speed (44 points of measure by 
rpm), the dispersive character of the central axes beam, 
compared to the stable mode, can be observed, where 
this same beam is tightened more and less tilted com- 
pared to the normal axis on the plane (x,y). This central 
axes dispersion can be explained by the self-excited vi- 
brations which cause the variable moment generation. 
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4.2 Analysis of central axis moments related 

While transporting the moment from tool tip to the cen- 
tral axis, the minimum moment (pure torque) Ma is 
obtained. From the moment values to the central axis, 
the constant and variable part of this one is deduced. As 
for the efforts, the variable part is due to the self-excited 
vibrations as revealed below (Fig. [T6|) . 

Using this decomposition, the moments contribution 
on the zones of contact tool / workpiece / chip is ex- 
pressed. The observations resulting from the analysis 
show that the tool vibrations generate rotations, cause 
variations of contact and thus generate variable moments, 
confirming the efforts analysis detailed in the section [27T1 
This representation enables to express the moments along 
the three axes of the machine tool: swivel moment in the 
y direction and the two moment of rotation along x and 
z directions. 



[clNin] 




[dNin] 



Fig. 16 Moment representation related of central axes; case 
study using ap = 5 mm, f=0,l mm/rev and N = 690 rpm. 

Moments components analysis determined at the cen- 
tral axis allows noting a moments localization mainly on 
two distinct zones. Taking into account this aspect, the 
variable moments components are divided into two noted 
parts di and d2 on to the three directions related to the 
machine. The components of these variable moments on 
x, y, z, axes are noted respectively (M aX v), (M ayv ) and 

(Mazv). 

In the vibratory case (ap = 5 mm), the first points 
family is located (Fig. [16]) according to a line cfo and 
d\ (the large ellipse axis). In the case without vibra- 
tions (ap = 2 mm) , the two d\ and ^2 families merged to 
only one family located according to only one line. Thus, 
the elliptic form appearance around the right-hand side 
d\ seems quite related to the self-excited vibrations (cas 
ap = 5 mm) [7]. In particular, the frequency associated 
with the d\ part is higher than those associated the ^2 
part. Furthermore, these frequencies are related to the 
frequencies domain found during moments FFT analy- 



sis Fig. H21 Finally, on the central axis, the d\ families 
and c?2 seem to correspond to distinct elements from the 
generated surface. 



5 Workpiece and chip geometry 

5.1 Roughness measurements 

In the processes of matter per cutting tool removal, it is 
well-known that the manufactured pieces surface quality 
is closely connected to the thrust force imposed on the 
matter [TT]. In particular, the more the thrust force ex- 
erted on the object surface is important, the more the 
imprint left by the tool on surface is consequent. The 
surface roughness east thus closely related to the thrust 
force intensity. The surface roughness is connected con- 
fidentially to the thrust force intensity. The self-excited 
vibrations have an influence on the workpieces quality 
surface. Now, in the self-excited vibrations case exam- 
ined here, we have just noticed that the forces and mo- 
ments are maximum with the resonance frequency which 
is situated around 190 Hz. 

Taking into account the link between roughness and 
efforts applied to the surface, we should find a maximum 
of surface roughness around this frequency. In this pur- 
pose, we propose to control the profile of a generator 
length roughness of the cylindrical manufactured piece. 
Indeed, in every rotation of the machined piece, the tool 
leaves a imprint force function applied to the generator 
point observed which corresponds to a given and known 
moment. Along a given generator, each point of this one 
is the force image applied to the surface at a known given 
moment because the rotation speed of the workpiece is 
known. 

The surface roughness examination along a generator 
should thus necessarily reveal the amplitudes of rough- 
ness associated with the efforts periodically applied by 
the tool in these points. Also the roughness profile FFT 
presents along a generator of the cylindrical part man- 
ufactured should pass, like the efforts, by a maximum 
around the frequency of 190 Hz. It is indeed what can be 
observed on the figure [TBI where the roughness data FFT 
analysis shows a frequency peak located around 190 Hz 
(precisely 191,8 Hz), which is coherent with the previous 
data. 

In addition, the surface roughness analysis gives a 
total roughness value R t = 1,6 fim. 

5.2 Chip characteristics 

Chip measurements under the Scanning Electron Micro- 
scop were carried out and enabled to determine the thick- 
ness variation and the chip width. All chips are type 1.3 
(ISO 3685) with undulations. 

The chip thickness variations between the maximum 
(h max ) and the minimal (hmin) thickness are about 2, 
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190 200 210 220 

Frequency 



Fig. 17 Machined surface and FFT data related to the 
roughness; case study ap = 5 mm, f = 0,1 mm/rev, N = 
690 rpm. 



and feed rate values independent. An example is pre- 
sented (figure H5j) for a sample of chip during a test with 
the feed rate f=0.05 mm/rev. Values obtained h max = 
0,23 mm and h m i n = 0,12 mm. 




Fig. 19 Chip width evaluation. 



[T9] of the chip during cutting process and the primary 
shear angle </>, 



cos 



(0-7) 



sin(0) 



(4) 



In our case, it is measured on the chip undulation 
length Iq = 11 mm, with a chip hardening rate £ c = 1.8 
and the cutting speed V = 238 m/min. A frequency of 
206 Hz is then obtained, very near to the frequencies of 
tool tip displacements or of the load application points 
during cutting process. 

The chip width is then measured with similar tech- 
niques. Substantial width variations are observed, about 
0.5 mm. Indeed, the measured maximum width Wm a x 
is 5.4 mm while the evaluated minimal width w m i n is 
4.9 mm, which means a width variation about 10% - 
Fig. [191 

The chip width slopes angle ip c measurement (Fig.fTO]) 
between each undulation, close to 29°, is equal, with the 
same errors, to the phase difference measured on the sig- 
nals of tool tip displacements (28°) (cf. [2]). 




Fig. 18 Chip thickness variation evaluation. 



6 Correlation between displacements of the tool 
tip / applied forces 

A synthesis between the work-paper two parts is essen- 
tial. It is carried out below in order to put in evidence 
the various correlations which exist between stiffness / 
displacements, displacements of tool / stiffness center or 
stiffness center / central axis. 



The measure of chip length corresponding to an un- 
dulation enables to find the self-excited vibrations fre- 
quencies starting from the cutting speed (in conformity 
with equation 3): 

V 

fcop = -j ^ 5 (3) 

with f cop the chip segmentation frequency, V chip 
speed, Iq one chip undulation length and £ c the chip 
hardening coefficient. 

To determine the total chip length, it is necessary to 
measure the wavelength, taking into account the rate of 
hardening phenomenon (in conformity with equation 4) 



6.1 Correlation between the plane of the displacements 
tool tip / the applied forces 

The tool tip point displacements are localized (cf. [4]) 
in a tilted plane. Inside the stiffness matrix determined 
in the paper (|6j), a correlation exists between tool tip 
displacements and the cutting forces applied. In partic- 
ular, the ratios between the large and the small ellipses 
axes of tool displacements (a u /b u ) and the efforts ap- 
plied (cLf/bf) decrease feed rate functions while the ratio 
of these ratios remains constant (equal to 1,64) when 
the feed rate increase. These elements enable to deter- 
mine accurately the real configuration of cutting process. 
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f (mm/rev) 


0.05 


0.0625 


0.075 


0.1 


Normal 


n f 


n u 


n f 


n u 


n f 


n u 


n f 


n u 


Along x 


0.245 


-0.071 


0.292 


-0.071 


0.419 


-0.058 


0.46 


-0.056 


Along y 


-0.107 


-0.186 


-0.113 


-0.186 


-0.097 


-0.206 


-0.1 


-0.216 


Along z 


-0.964 


0.98 


-0.95 


0.98 


-0.903 


0.997 


-0.882 


0.975 



Table 4 The normals (nf ,n u ) of the tool points displacements planes and the forces applied along the three cutting directions, 
depending on feed rate. 



These correlations are analyzed using the direct normal 
to the tool tip point displacements plane and the direct 
normal to the load application points place (Table [4)) . 

The existence of these two planes is particularly inter- 
esting and adapted to the establishment of a cutting pro- 
cess real configuration. This new aspect is in the course 
of implementation in order to express and exploit a sim- 
plified dynamic three-dimensional model in the reference 
system associated with these planes [5]. 



6.2 Self-sustained vibrations: experimental validation 

From these studies, it comes out that the self-excited vi- 
brations domain is around 190 Hz with accuracy of a 
few per cent. It is around this common fundamental fre- 
quency that the whole of the major characteristics (dis- 
placements, efforts, and moments) of system WTM have 
the most important amplitude variations. The analysis 
carried out to the measures of tool point displacements 
and the points for load application enables to evalu- 
ate existing constant phase difference between the forces 
components and corresponding displacements (Fig.l20lto 




Time (a) 



Fig. 20 Phase difference evaluation between forces / dis- 
placements along x axis. 



These phase differences confirm the efforts delay com- 
pared to the tool tip displacement. The self-excited vi- 
brations appearance can be also explained by the delay 
force / displacement, which increases the system energy 
level. The existence of this delay could be explained by 
the machining system inertia and more particularly by 
the cutting process inertia [18J. 
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Fig. 21 Phase difference evaluation between forces / dis- 
placements along y axis. 




Time (s) 



Fig. 22 Phase difference evaluation between forces / dis- 
placements along z axis. 









13° 


23" 


75" 



Table 5 Values of phase difference forces / displacements 
(ap = 5 mm, f = 0,1 mm/rev, N = 690 rpm.) 



Moreover, it is noted that the phase difference with 
the same range between efforts / displacements (Table EJ) 
remains constant according to the feed rate. 

Because of parts elasticity intervening in the opera- 
tion of turning, it is logical that the response in displace- 
ments of unit BT-BW is carried out with a certain shift 
compared to the efforts variation applied to the tool, 
variation induced by the lacks of machined surface cir- 
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cularity in the preceding turn which imply variations of 
the contact tool/workpiece (Fig. [17)). Phase difference 
between the efforts and displacements thus remains a 
possible explanation to the self-excited vibrations ap- 
pearance. 




Fig. 23 Movement tool / workpiece elliptical trajectory. 

Moreover, when the tool moves along the ellipse places 
ei, 62, es (Fig. [23]), the cutting force carries out a pos- 
itive work because its direction coincides with the cut- 
ting direction. On the other hand, on the side es, e^, 
ei, the work produced by the cutting force is negative 
since its direction is directly opposed to that of displace- 
ment. The comparison between these two ellipse parts, 
as shows that the effort on the trajectory ei, e2, es is 
higher as on the trajectory es, e±, e\ because the cutting 
depth is more important. At the time of this process, 
work corresponding to an ellipse trajectory remains pos- 
itive and the increase of result energy thus contributes 
to maintain the vibrations and to dissipate the energy in 
the form of heat by the assembly tool/workpiece. 



a plane inclined compared to the machine tool reference 
system. 

This cutting force, whose application point describes 
an ellipse, is perfectly well correlated with the tip tool 
displacement which takes place under similar conditions. 
In particular, the ellipses axes ratios remains in a con- 
stant ratio when the feed rate varies with a proportion- 
ality factor equal with 1,64. 

Moreover, the highlighted coupling between the sys- 
tem elastic characteristics BT and the vibrations gen- 
erated by cutting enabled to demonstrate that the self- 
excited vibrations appearance is strongly influenced, as 
it was expect to, by the system stiffness, their ratio and 
their direction. We also established a correlation between 
the vibratory movement direction of the machine tool 
elastic structure, the thickness variations and the chip 
section. 

These results enable to now consider a more complete 
study by completely exploiting the concept of torsor. In- 
deed, thanks to the six components dynamometer, we 
confirmed, for an turning operation, the moments ex- 
istence to the tool tip not evaluated until now by the 
traditional measuring equipment. 

The originality of this work is multiple and in partic- 
ular consists in a first mechanical actions torsor analysis 
applied to the tool tip, with an aim of making evolve 
a model semi analytical cutting 3D. This study allows, 
considering a turning operation, to establish strong cor- 
relations between the self-excited vibrations and the me- 
chanical actions torsor central axis. It is thus possible, 
thanks to the parameters use defining the central axis, 
to study the vibrating system tool/workpiece evolution. 
It also leads to the description of a "plane of tool tip 
displacements" in correspondence with "the load appli- 
cation points plane". 

Thus, using the plane that characterizes the BT be- 
haviour enable to bring back the three dimensional cut- 
ting problem, with space displacements, with a simpler 
model situated in an inclined plane compared with the 
reference system of machine tool. Nevertheless, that re- 
mains a specific model of three dimensional cutting. 



7 Conclusion 

The experimental procedures proposed in this work-paper, 
as well at the static and dynamic level, enabled to de- 
termine the elements necessary to a rigorous analysis of 
the tool geometry influence, its displacement and evolu- 
tion of the contacts tool/workpiece and tool/chip on the 
machined surface. 

In particular, by analyzing the efforts resultant tor- 
sor applied during turning process, the experimental re- 
sults enabled to establish an efforts vector decomposi- 
tion highlighting the evolution of a variable cutting force 
around a constant value. This variable effort evolves into 



8 Appendix 

8.1 Determination of the place points plane of load 
application on the tool 

The plane P u definition being load application points ge- 
ometrical place on the tool starting from the experimen- 
tal results is carried out using the computation Mathcad© 
software. We seek to determine the plane which passes 
by the load application points cloud on the tool (Fig. [8] 
section [2^4|) : 

ax + by + cz + d = 0. 
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The errors are noted with e rr and we have : 



6 rr (x, UiZ^ l/p, Zp, X ni y<m Zn) — 

[M(x, y, z) - P(x p , y p , z p ).n f (x n ,y n , z n )] , 
where: 

M = {x, y, zf , P {x p , y p , z P ¥ , n f = {x n , y n , z n } 1 . (5) 

Here the superscript t indicates the operation of trans- 
position. 

Expressing the E rr function using e rr and introduc- 
ing the displacements components (u Xi ,u yi ,u Zi ) of the 
tool load application points into the three space direc- 
tions, it comes: 



E rr {x^ y z , x p: yp-, Zp^ x n: ynt Zn) — 
^ ^ Crri^Xi •) Uy i , u Zi , Xp, yp, Zp, x ni y ni z n ), 



TV 



i=0 

and : 

EN sr^N sr^N 

i=0 U Xi _ l^i=Q u Vi _ 2-^i=Q u Zj 

Xp ~ N " Vv ~ N ' p ~ N ' 

Now, the vector V is calculated by minimization: 

V = TniTi(E rr , Xp : yp^ Zp, x n: ym Zn)i 
where 



Pf(s, t) = V P + s.ui + t.u 2: 

where s and t are constants, and V p is equal to 




Ui 



and ui is given by: 
uo A nf 

= Ti a N' U2 = 11 AUl ' 

||n A u f || 

with the plane uq orientation vector is: 



u 



8.2 Ellipse approximation 

Using the ellipse plane determination [3j, it is possible 
to determine the characteristics given in Table 03 
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r v 1 




v 2 


= y P 


< v * 


%p- 


v 4 


x n 


v 5 


= y n 


{ V 6 


— z n 



It results from it that the direct normal rtf compo- 
nents to the plane of the load application points on the 
tool are given by: 

f V 4 
n f = ^ V b . 

{ v 6 

For the case study presented in the section 12.41 (Ta- 
ble [2]), with ap = 5 mm and f= 0.1 mm it comes: 



+0.46 
L f = <; -0.1. 

-0.882 

The required plane equation is then: 
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